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ABSTRACT 


Analysis has been performed for mapping VISSR imagery 
data so as to eliminate all geometrical distortions. The formula- 
tion is rigorous and includes all misalignment angles of the VISSR, 
the sun sensor and the instantaneous spin axis with the satellite's 
body axis. It also includes the effects due to the motion of the 
satellite's suborbital point. All the mapping equations for dis- 
tortion removal are reduced 'to simplest forms, and all the algorithms 
are optimized as much as possible. 

An approach is then formulated for implementing these 
algorithms for in-line operational use . It covers the computations 
involved in determining benchmarks, the interpolation methodology 
for filling in the points interspaced between benchmarks, and the 
correction procedure for computing the radiranetric values at the 
center of the pixel in the distortion-free image. It is also concerned 
with the time requirements, data storage, and output data accuracy. 
With the present microprocessor technology, it is concluded that this 
in-line distortion removaJ. is possible in real-time processing of 
infra-red but not visible VISSR imagery data* 


This work was supported by NOAA Contract Nos; OI- 8 -MOI-I 86 U 
and NA-79KAC-00026 
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SECTION 1 - INTRODUCTION 

In the Visible and Infra-red Spin Scan Radiometer (VISSR) 
data obtained from the present geosynchronous satellites, distortions 
are obsei*ved in the images of the earth. As illustrated in Figure 1 .1 
which is exaggerated for clarity, these image deformations appear 
as vertical compression and expansion of the image, non-vertical 
alignment of the North and South Poles, and multi-representation of 
some points or omission of other points. 



Figure 1.1 Exaggerated VISSR Image of Earth 
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The3d distortions may well be explained by considering 
Figure 1,2 which illustrates the same scan-lines on the projection 
plane of the earth. Again, for clarity, these scan-lines are depicted 
to be non-parallel and unevenly spaced to a degree more so than the 
realistic cases. 



s 


Figinre 1,2 Projection Plane Image of Earth 
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If one were to relate these two images, one would find, 
for example, that a triangular figure in the projection plane 
image becomes distorted into a curved figure in the VISSR image* 

This is illustrated in Figiire 1,3 which is obtained by super- 
imposing a triangle on Figure 1.2 and then mapping it onto Figure 1.1, 
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Projection Plane Image 




VISSR Image 


Figure 1,3 Distortion of Image 
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The gener&l causes for these image deformations may be 
broadly classified as follows: 

1. Orbit not circular and equatorial 

2. Spin axis not perpendicular to orbital plane 

3. Misalignment of the VISSR, the sun sensor, and the 

instantaneous spin axis with the satellite's body axis 

li« Biases due to varying sun size and varying s\in elevation 
effect on threshold of the siai sensor triggering. 

To remove these distortions, it is necessary to include all the 
above factors in the formulation of the mapping equations. However, 
it is feasible to consider only the first three. The corresponding 
equations have been derived in Reference 1 in which it was convenient 
to introduce the following coordinate systems : 

The Inertial System : This is well-known and is defined such that 

the Xj-axis is in the direction of the vernal equinox, the Zj-axis 
is perpendicular to the equatorial plane (in the direction of the 
North Pole), and the y^-axis is given by y^ » z"j x x^. 

The Body System : This system is defined such that the Zg-axis is 

along the longitudinal axis of the satellite, the Xg-axis is the 
intersection line between the VISSR stepping plane and the plane 
perpendicular to the Zg-axis, and the yg-axis is given by yg ■ Zg x ^ 
The VISSR System ; This system is defined such that the x^-axls is 
in the direction of the mid-scan, the z^-axis is perpendicular to the 
x^-axis and lies in the VISSR stepping plane, and the y^-axis is 

. . A /S 

given by y^ ■ z^ x x^. 

The Sun Sensor System : This system is defined such that the Xy-axis 
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is in the direction of the bisector of the anglh fanned out by the 


sun sensor, the z^-axis is perpendicular to th$ x^-axis and lies in 
the sun sensor fan plane, and the y^-axis is given by y^ ■ 'zy x 
The Spin System ; Let Zg denote the \mit spin axis vector around 
which the satellite is instantaneously rotating. Let r' denote the 


position vector of the satellite. Then, the y_-axis and the Xg-axis 


are respectively defined by 



A 

% = 


A A 

% ^ h 
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The Aiixiliary System ; In this system, illustrated in Figure l.li, 
the \init base vectors are defined by the following equations : 
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The Normalized System ; Consider a system, referred to as the normalized 
system, as illustrated in Figure f,$. The origin R of this system is 
defined at a point"? on the earth’s equatorial plane and fixed in the 
earth’s rotating system. The satellite P at point however, is not 
necessarily on the earth's equatorial plane or fixed in the earth's 
rotating system. In this normalized system, the unit base vectors 


are defined by the following equations: 


A A 

ss "A 



A A 
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In the present study, the res^llts obtained in Reference 1 
are used to formulate algorithms for mapping the data coming out of 
the Synchronous Data Buffer (SDB) so as to obtain a distortion-free 
imaigery. Moreover, this rectified imagery eilso has the desirable 
feature that it is referred to a normalized satellite position which 
is therefore the same for all imageries. Thus, if the distortion-free 
mapping is performed in-line during data processing, the transmitted 
VISSR data will provide a uniformly compatible data base for all 
users in their scientific work. Furthermore, it will also facilitate 
in the future development of a composite data base for different 
kinds of data obtained frcm various satellites. 

Mapping of the data may be further optimized the use of 
interpolation with the aid of appropriately chosen benchmarks. 

Section 2 deals with the computation of these benchmarks, while 
Section 3 covers the interpolation methodology for filling in the 
points interspaced between benchmarks . Section U is concerned with 
the correction procedure for computing the radiometric values at the 
center of the pixel in the distortion-free image. Section 3 discusses 
the time requirements, data storage and output data accuracy. 

Section 6 summarizes the results of this study. 
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SECTION 2 - BENCHMARK COMPUTATIONS 


This section deals with the computation of benchmark 
locations in the normalized distortion-free coordinate system. It 
discusses the relevant input parameters, computational equations, 
number of corapuational operations, and the requisite partial deriva- 
tives. 

2.1 Input Parameters 

The relevant input parameters are listed below ; 

N^ = scan number 

N, = sample number 

= mid-scan, number 

Mp = mid-sample number 
>s 

A, = scan angialar width 
= sample angular width 
“ scan angular bias (line bias) 
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"r* ® orbital position of Satellite 


a - right ascension of spin axis 


S « declination of spin axis 


"r = normalized position of satellite 


a = earth's semi -major axis 
c = earth's semi -minor axis 

For convenience, two parameters dependent on tne above are defined 
as follows: 
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2.2 Computational Equations 


The following equations for computing benchmark coordinates 
have been extracted from Reference 1 . They have been simplified and 
are listed below in 'the proper sequence for usage. The exact definitions 
of cursory intermediate variables may be obtained from the original 
report. 
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j* = 


1“ j’"' are the coordinates in the normalized distortion-free 

system. At this stage, for the sake of greater accuracy in sub- 
sequent computations, it is preferable not to digitize them. 
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2.3 Computational Operations 

For each benchmark, it may be verified that the comput- 
ations in equations (2.1 ) - (2.28) may be achieved by performing: 

additions 
22 divisions 
1 1 ii multiplications 
ii divisions 

1 1 trigonometric function evaluations 
2 square root evaluations 

Assuming that the following times are required: 


Operation 

Time (micro- 

Addition 


Subtraction 

1.5 

Multiplication 

6 

Division 

11 

Trigonometric function evaluation 

50 

Square root evaluation 

50 


it is seen that about 1U87 microseconds are required for each bench- 
mark computation. 
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2.h Partial Derivatives 


For convenience, it is desirable to choose the set of 
3 so that they form a rectai 
-space as illustrated in Figure 2.1. 


benchmarks so that they form a rectangular grid in the (n^ , N ) 



Figure 2.1 - Benchmarks in (N_ > N ) - Space 

X 


"X** •){■ 

Then, it Is obvious that the partial derivatives of I and J 
with respect to N> and N may be approximated by 


•7 



(i-j - (i‘), 
(«,. ), - ), 
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ji")^ - (j% 
(r)n, - (r)p 


(2'3f) 




For each benchjnark, the four associated partial derivatives 
require 6 subtractions and h divisions. These operations consume 
about 53 microseconds. 


These partial derivatives are used later in the method of 
interpolation for mapping points interspaced bet^^eeen the benchmarks. 
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SECTION 3 - INTERPOLATION COMPUTATIONS 

This section is concerned vdth the mapping of points 
which do not coincide with the chosen set of benchmarks. It 
discusses the input data, interpolation methodology, and number of 
computational operations. 

3.1 Input Data 

The input data consists of the coordinates (T^, and 
their four associated partial derivatives for each benchmark. This 
information has already been obtained in Section 2. 

3.2 Interpolation Methodology 

Suppose there are m interspaced points between the hori- 
zontal benchmarks, and n interspaced points between the vertical 
benchmarks. Figure 3«1 illustrates a basic unit comprising bench- 
marks (denoted by solid circles) and interspaced points (denoted by 
open circles). Thus, there are (m + 1)(n +1) points altogether 
in a basic unit. 
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Figure 3*^ " Benchmarks for Interpolation 


It is noted that the I" coordinate of the point T directly 
below P is given by 
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Therefore, it follows that 


(r)r = + l(^)J 


(a- 4) 


, 'Wf/ P 


A similar equation holds for any point and the point directly 
below it. Hence, the n interspaced points between P and R may be 
mapped by the following iterative algorithm; 


Let 




C3>s-y 


j: - o*\ 


Ci'i) 


Then, perform the following computations 


I 


Ml 


J 


C-tl 


r«- 


T* 


[(S^ ij 


f/iia 


,u 




0,1,2^. ..n 




Similarly, it is noted that the I” coordinate of the 
point U directly to the right of P is given by 


(z^ = 
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Therefore, it follows that 




r4 - (i% + 

Again, a similar equation holds for any point and the point directly 
to the right of it. Hence, the m(n +1) interspaced points between 
the columns PR and QS may be mappea by the following iterative 
algorithm: 


Let 



( 3 >/ 3 ') 

o, 1,1^ ... n 


Then, perform the following computations 
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After these computations, it is now permissible to digitize I and J 
to obtain the rounded integer values I and J. 
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3.3 Computation Operations 


For each interspaced point, it is seen that the comput- 
ations in equations (3.7) and (3.8) or those in equations (3.15) 
and (3.16) require 2 ’’additions. Assuming a time of 1.5 microse- 
conds for each operation^ -therefore about 3 microseconds are required 
to map each point by interpolation. Allowance is also to be make for 
converting tv;o real numbers to integer values for each point. This 
will probably increase the time requirement by a factor of 2 so 
that about 6 microseconds are required for each point. 
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SECTION li ~ RADIOMETRIC COMPUTATIONS 

This sections considers the methodology of correcting the 
radiometric values so as to reflect a more realistic value at the 
center of the pixel in the distortion-free image. It discusses the 
input data, correctidn methodology, and number of computational 
operations. 

it.1 Input Data 

For each point, the input data consists of the follov7ing 
and may be obtained either from the SDB output data stream or has 
already been obtained in Section 3: 


N„ = scan number 

’I 


N- = sample number 

h 


R(N^ , ) = radiometric value of th scan and th 

sample 



“ radiometric value of th scan and 

(N - - 1 )th sample 

li 


R(N j ”1 ) = raidometric value of (N^ - 1 )th scan 

and Nj, th sample 

is 
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I" * interpolated horizontal coordinate, -of (N^ , ) 

in distortion-free system 


J" “ interpolated vertical coordinate of (Nw. , N ) 
in dis’tortion-free system 


I = rovinded integer value of 1' 


J = rounded integer value of 


U.2 Correction Methodology 


The radiometric value R(I, J) may be obtained from the 
value R(I", J") by using the Taylor's series expansion 

ft X) - X') X fi-i') t 

(4'l) 

f \ fL\ 

The partial derivatives f J and "i^y be written as 

T* I* 






r 


^ fiL ) (3i] + (2L) m 


J* 


) ^ /'2£-) ^ flS.) 

x' WxV;j. 




(4S) 



The partial derivatives (_ 9 /^ i ’and j "lay be approxi- 
mated by ^ % 









The four remaining partial derivatives on the RH3 of equations [h»2) 
and (U.3) may be obtained as follows: Let A emd B be matrices defined 

by 
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Then, from the theory of mathematicaX transformations, we have 

h = fi”' 


which explicitly yields 





I ni*] 
C ( 




where 



The partial derivatives appearing on the RHS of equations (U.9) - 
may be obtained from equations (2.29) - 12.32), valid for ; 
basic unit defined by benchmarks. 
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U. 3 Computational Operations 


Jj’or each point, it is seen that the computations in 
equations (U. i / - (U»5) require I4 additions, 6 subtractions and 
6 multiplications. Assuming a time of 1.5 microseconds for each 
addition or subtraction , and 6 microseconds for each multiplication, 
therefore about 51 microseconds are required to correct the radio- 
metric value for each point. 
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SECTION ^ - DISCUSSION 

This section discusses relevant topics such as benchmark 
spacing, time requirements for mapping IR data, real and non-real time 
computations, data storage and buffering, input and output data ac- 
curacy, and computational accuracy requirements* 

5.1 Benchmark Spacing 

The IR samples have angular widths of about 0.01° x 0.005° 
at the satellite position. This corresponds to a resolution of about 
U X 2 miles at the sub-satellite position on the earth's surface. In 
general, this resolution and the non-linearity of the mapping equations 
determine the requisite spacing of the benchmarks to be used for inter- 
polation. The analytical approach to obtain this spacing involves 
comparatively complex mathematical anailysis. Alternatively, it is 
also possible to obtain this value by actually performing the mapping 
munerically. At this stage, it is felt that the interpolation require- 
ments can be met by choosing the IR benchmarks to be spaced 50 samples 
horizontally and 25 samples vertically. That is, it probably suffices 
to choose m = 50 and n » 25 in Section 3. In the full IR imagery, 
there are 1322 scans each containing 3822 samples. Consequently, about 
5,600 benchmarks will be required. 
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5.2 Time Requirements 

In Section 2, the computation of each benchmark and its 
associated partial derivatives requires about ^ShO microseconds. 

Hence, a set of $S0Q benchmarks requires about 8,62h,000 microseconds 
— 8.6 seconds. 

In Section 3> the mapping of each sample by interpolation 
requires about 6 microseconds. Therefore, an IR imagery of about 
7 X 10^ samples requires about 1|2 seconds. However, if the entire 
IR imagery is not to be mapped, then cropping out the edges will 
probably reduce time by a factor of 2/3 to yield a requirement of about 
28 seconds. 


In Section U, the correction of radiometric value at 

the center of the pixel in the distortion-free image requires about 

microseconds for each sample. Therefore, an IR imagery of about 

6 

7x10 samples requires about 357 seconds. Cropping will probably 
reduce this to about 238 seconds. 

Consequently, about 399 seconds or 6.7 minutes will be 
needed to map the entire IR imagery crsmprising of coordinates and 
radiometric values of the samples. This time requirement drops to 
about U.U minutes if cropping is introduced. 
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If it is desired to map»the visible imagery containing 
1/2 X 1/2 mile samples, then the above times are increased by a 
factor of 32, Therefore, about 21 U minutes will be required to map 
the entire imagery comprising of coordinates and radiometric values. 
If the edges are cropped out, then about 1U3 minutes will be needed. 

^.3 Real and Non-Real Time Computations 

From the discussion above, it is seen that it is possible 
to perform all the mapping computations in real-time in the case of 
IR imagery, and not possible in the case of visible imagery. However 
in the latter case, the crucial point is whether the radiometric 
corrections are really necessary. If not, then the time requirements 
drops to 22. U minutes for the entire imagery, and lU.9 minutes for 
the cropped imagery. Consequently, visible data-mapping becomes feas 
ible in real-time. 

Because the benchmark computation time is so small, it 
is desirable to perform the benchmark computations in real-time so 
that the relevant parameters may be easily extracted in-line from 
the data-stream coming out of the SDB, 
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$.h Bata Storage and Buffering 

Because the imagery obtained from the SDB output data is 
distorted, it is necessary to store this data in a buffer before the 
distortion-free mapping can be performed. The buffer size may be 
estimated by allox-ring for a maximum 3° offset in the spin axis. Since 
the satellite is about 6.6 earth radii away, it may be verified that 
about 100 IR scan-lines (382,200 samples) to be buffered at a time. 
This wfill be sufficient to output a horizontal distortion-free line 
from end to end. In the case of visible data, the corresponding 
buffer will contain about 800 visible scan-lines (12,230,[t00 samples). 
If a realistic situation, the above numbers will probably be reduced 
by a factor of 3» 

5.5 Input and Output Data Accuracy 

The data coming out of the SDB will be used as input 
into the distortion-free mapping software system* The acciiracy of 
this data may be roughly classified as perfect, normal or bad. 

Perfect data corresponds to data having errors of less 
than one pixel (i. e. ♦ 2 km at the subsatellite point for IR data). 

The error in the output data from the distortion-free mapping is 
therefore determined by the pixel resolution of the benchmarks, the 
interpolation accuracy/ of interspaced points, and the correction accur- 
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acy of the radiometric values. Because the second and third factors 
depend on the benchmark-spacing, which in turn de^)ends on the pixel 
resolution, therefore it is estimated that the error bound of the 
output data is about one pixel (i.e., + U km for IR data). 

Normal data corresponds to data having errors of about 
one or two pixels. The mapping error is determined by the benchmark 
accuracy corresponding to normal input error, the interpolation 
accuracy of interspaced points, and the correction accuracy of the 
radiometric values. In this case, the error bound of the output 
data is about two pixels. 

Bad data corresponds to data having errors of about U or 
more pixels. The mapping error is determined mainly by the benchmark 
accuracy corresponding to these bad input errors. In this case, the 
error of the output data is probably about ^ or more pixels. 

5.6 Computational Accuracy Requirements 

It is desirable to investigate into the use of 1 6-bit 
words in the distortion-free computations. 

Because of the complexity of the benchmark computations 
in equations (2.1) - (2.28), it is quite evident that sufficient 
accuracy will not be obtained by performing single-precision 
computations using 1 6-bit words. 
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However, for the partial derivatives comnutations in 
equations (2.29J - (2,32;, the interpolation of interspaced points 
computaions in equations (3.7; - (3.0) and (3.15) - (3.16), and the 
radiometric correction computations in equations (U.i; - kU.5), 
it is possible i<o achieve the desired accuracy using single-precision 
computations involving 1 6-bit words. In this case, perhaps the best 
way to represent real numbers is as follows: 

I bit for sign of niimber 

1 1 

I I bits for range of niomber ( 2 

1 bit for sign of exponent 

3 

3 bits for range of exponent ( 2 
An alternative choice is as follows: 

1 bit for sign of number 

12 

12 bits for range of number (2 - 1 = U095 ) 

1 bit for sign of exponent 

2 

2 bits for range of exponent (2 -1=3) 

This second choice may not be as desirable because of the small 
range of the exponent. 


- 1 = 20U7 ) 
-1=7) 
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SECTION 6 - CONCLUSION 


PVom the preceding discussion, it is seen that it is 
possible to map in real-time the entire IR imagery comprising of 
coordinates and radiometric values of the samples. However, it is 
possible to map in real-time the entire visible imagery comprising of 
only the coordinates of the samples. This conclusion is based heavily 
on the assumption that it takes 1.5 microseconds for each addition or 
subtraction, and 6 microseconds for each multiplication. 

The follox«.ng table summarizes the time requirements for 
IR and visible imagery mapping: 



Entire IR 

Cropped IR 

Entire VIS 

Cropped VIS 

Benchmarks 

8.6 sec. 

5.7 sec. 

U.6 min. 

3.1 min. 

Sanole 

Coordinates 

U2 sec. 

28 sec. 

22. U min. 

1h.9 min. 

Radiometric 

Values 

3^7 sec. 

238 sec. 

I9O.U min. 

126.9 min. 

Total 

6.8 min. 

U.5 min. 

217. U min. 

IUL.9 min. 
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The following table summarizes the expected accuracy of 
the distortion- free mapping (DFM) algorithms: 


Input Data 

Output Data 

from SDB 

from DFM 

Perfect 

1 pixel 

Normal 

2 pixels 

Bad 

5 or more pixels 
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